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Based on the most recent data in charmless B decays including the very recently reported large 
direct CP violations, it is shown that the weak phase 7 can well be extracted without two-fold 
ambiguity even only from two decay modes 7r"'"7r~ and 'k'^ K~ , and its value is remarkably consistent 
with the global standard model fit at a compatible accuracy. A fit to all the -k-e^-eK data favor 
both large electroweak penguin and color-suppressed tree amplitude with large strong phases. It 
is demonstrated that the inclusion of SU(?>) symmetry breaking effects of strong phases and the 
inelastic rescattering effects can well improve the consistency of the data, while both effects may 
not be sufficient to arrive at a small electroweak penguin amplitude in the standard model. It is of 
interest to notice that large or small electroweak penguin amplitude becomes a testable prediction 
as they lead to significantly different predictions for the direct CP violations for tt^tt", if'K'^ modes. 
Clearly, precise measurements on charmless B decays will provide a window for probing new physics. 

PACS numbers; 13.25.Hw, ll.30.Er, ll.30.Hv 



The evidences of direct CP violation in B decays have 
recently been reported by the BaBar and Belle collab- 
orations. The latest averaged data give acp (7r"'"7r~) = 
0.46 ± 0.13 and acp(7r+A'-) = -0.11 ± 0.02[i|. Thus 
direct CP violation has been established not only in the 
kaon system, but also in B system. It has been shown 
that both direct CP violation and A/ =1/2 rule in kaon 
decays can be understood in the standard model (SM)|3]. 
It is then natural to test whether the observed direct CP 
violations and decay rates in charmless B decays can be 
explained within the SM. As the two experimental groups 
BaBar and Belle have reported more and more accurate 
data for charmless B decays {B it then 

allows one to test the SM and to explore possible indi- 
cations for new physics, such as the two-Higgs-doublct 
model with spontaneous CP violation^, the supersym- 
mctric models etc. There have been several global ana- 
lyzes which are based on either model independent pa- 
rameterizations such as quark flavor diagrammatic de- 
composition 0, 1^ QJ3i isospin decomposition^, fla- 
vor SU(3) symmetry [llj. or QCD inspired calculations 
such as QCp factorization and perturbative QCD 

approach[l3| as well as soft-collinear effective theory [lj|. 

In this letter, we are going to make a step-by-step fit 
for the charmless B decay modes based on approximate 
SU(3) flavor symmetry and hierarchical structures of di- 
agrammatic amplitudes. Based on the most recent data 
including the very recently reported large direct CP vio- 
lations, we arrive at the following main observations: i) 
the current data allow us to precisely extract the weak 
phase 7 from only two modes tt+tt^ and -n^ K~ without 
two-fold ambiguity. The resulting numerical value of 7 
is found to be remarkably consistent with the global SM 
fit at a compatible accuracy; ii) A direct fit to all tttt. 



ttA' modes favors a large electroweak penguin. Further- 
more, the large or small electroweak penguin amplitude 
is found to be a testable prediction via measuring direct 
CP violations in the decay modes B Tr°n° and tt°K° 
once more accurate data become available, iii) all the 
amplitudes and strong phases in _B ^ tttt, ttK arc ex- 
tracted, which indicates large final state interactions and 
non-factorizablc QCD effects as the resulting numerical 
results show an enhanced color-suppressed tree ampli- 
tude and strong phase. It is shown that not only the 
large tt^tt branching ratios but also the large tt^X^ ones 
result in a large color-suppressed tree amplitude with a 
large strong phase; iv) it is the large n'^K branching ra- 
tio that mainly responsible for a large electroweak pen- 
guin amplitude with a large strong phase. In the case 
of a small electroweak penguin amplitude fixed by the 
isospin relation in the SM, the resulting ir'^K branching 
ratios arc bellow the experimental data; vi) SU(3) sym- 
metry breaking of strong phases and B DD rescatter- 
ing effects can well improve the consistency of the global 
fitting. However, it remains necessary to have a large 
electroweak penguin amplitude with large strong phases. 

The diagrammatic decomposition approach is adopted 
to carry out a global analysis. The advantage is that 
in such an approach some decay modes can form, in a 
good approximation, closed subsets, which allows us to 
determine the relevant parameters without knowing the 
others. Although the number of data points decrease for 
each subset, the number of the free parameters decrease 
as well. Qf interest, the precision of the determinations 
is not necessarily lower than that of the whole global fit. 
Furthermore it may avoid the complicity and the poten- 
tial inconsistency in the current data when more decay 
modes are involved in the whole global fit. The compar- 
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ison between different results from different subsets may 
provide us important hints to understand those decays. 
In general, all the B tttt decay modes can be written 
in terms of diagrammatic amplitudes: tree (T), color- 
suppressed tree (C), QCD penguin (V), electro- weak 
penguin (Vew)i color suppressed electroweak penguin 
( V^^,) etc. The corresponding diagrams in B ^ ttK 
are denoted by primed ones, such as T', T", etc. Using 



the CKM factors A, 



(^) 



V* 



■\q — yqd{s)^ib, and the unitarity of 
the CKM matrix, the penguin type amplitude can be de- 
composed as : T'C) = A^fV^ + xi'^'^Vc + >H^Vf Defining 
P = Ptc = 'Pt — 'Pc , Ptu ^'Pt — Pu, Pew = Pew + Pew 
and factorize out the CKM factors, we arrive at the fol- 
lowing diagrammatic decomposition 
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where the rescaled amplitudes have a hierarchical struc- 
ture T ^ P ^ Pew- The primed and unprimed 
amplitudes are equal in the SU(3) limit. For simplic- 
ity, throughout this paper, we will neglect the smallest 
amplitudes of Pew take in a good approximation 
PEW,tu — Psw.tc = Pew and Ptu — Ptc — P due to 
t-quark dominance. As a phase convention, we take T to 
be real, i.e. 6t = 0. The amplitudes arc normalized to 
the CP averaged branching ratio Br = {\A\'^ + \A\'^)/2 
in units of 10~^, where the tiny differences due to the 
B° and B^ lifetime difference and the final state phase 
spaces are neglected. The direct CP violation is de- 
fined through acp = - + \A\^). The 
flavor SU(3) symmetry breaking effects for amplitudes 
are considered as \T'/T\ = \P'/P\ = \Pew/Pew\ - 
Ik/ fir — 1-28 from naive factorization. The SU(3) sym- 
metry breaking effects of strong phases are characterized 
by the phase differences of the primed and the unprimed 
amplitudes AS a = — Sa with A denoting for any of 



the amplitudes T, P, Pew etc. 

The decay modes of tt'^tt" and tt^ provide five data 
points: two CP averaged branching ratios _Br(7r+7r^) = 
4.6 ± 0.4 and Br{TT+ K') = 18.2 ± 0.9, two direct CP 
asymmetries acpij^'^i^^) and acp{T^~^ K~), and one mix- 
ing induced CP asymmetry S.,^-^ = —0.61 ± 0.14. Taking 
the flavor SU{?>) relations and neglecting Pew^ ^'^^ 
decay modes only involves T, P, Sp and the weak phase 
7. Thus all of them can be directly determined. A fit to 
the current data gives the following results 



in = 0.53 ±0.03, 



Sp = -0.48: 



-0.09 



|P| = 0.09 ±0.002, 

- 1 1-1+0.11 

7 — -L- J- J- -0.14 



(2) 



with a Xmin/'^-'^-f = 0.71/1. Where the well measured 
result of sin 2/? 0.73 ± 0.037 from B ^ J/ipKs has 
been used to relate the weak phase a to the weak phase 
7 via unitarity relation. The values of \T\ and |P| are 
well determined with relative errors less than 10%. The 
error of (5p is larger but can be reduced with more accu- 
rate data in the recent future. The ratio |P/P| is found 
to be around 0.17. Note that the best fitted angle 7 is 
in a remarkable agreement with the one from the global 
SM fit of the unitarity triangle which gives 7 = l.OSlg ^i 
and at a compatible accuracy. We emphasize that the 
above results arc obtained without the interference with 
other TTTT, ttK modes in which more diagrammatic param- 
eters C and Pew are involved. Therefore it provides a 
very promising way to extract 7 from charmless B decays 
and an important reference point for any further analy- 
sis. In obtaining the above result, the newly reported 
acpiT^^ K^) plays a key role. Without it, as shown in 
ref[^ , the determination of 7 suffer from a two-fold 
ambiguity with the other solution at 7 ~ 40°. In Fig^ 
we plot the Xmin as a function of 7. It is clearly seen that 
after including acp(7r^A'^) the global minimum (best- 
fit) of falls into the allowed range of the SM fit and 
the ambiguity is lifted. 

Note that in the above fit the tttt and ttK modes are re- 
lated via the SU(3) relations, while the symmetry break- 
ing effects on the strong phases have been neglected. 
As pointed out in ref. |9( the SU(3) breaking in strong 
phases may significantly change the correlation between 
acp(7r^7r~) and acp(7r^i^~). Taking ASp as a free pa- 
rameter in the fit, we find 



\T\ = 0.53 ±0.03 

Sp = -o.67l°:^t 

7 1.06 ±0.2 



|P| = 0.09 ±0.002 
ASp ^ 0.21 ±0.4, 



(3) 



with xLn = 4.7 X 10 



which manifests that a small 
value of ASp further improves the goodness-of-fit. 

When including the branching ratios of Pr(7r°7r°) = 
1.51 ±0.28 and Br{-K°TT~) = 5.5 ±0.6 but ignoring P^vk 
at the moment as both modes are dominated by T and 
C, only two new parameters C and Sc are involved. A 
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FIG. 1: The Xmin as functions of 7. The three curves (from 
bottom to top) are: SoUd: Fit to tv^tt^ and tt'^K^ data only. 
Dashed: Fit to all the tvtt, tvK modes, with Pew free (Fit 
B). Dotted: Fit to all the nn, ttK modes, with Pew fixed 
by Eq.(|5Jl (Fit A). The shadowed band indicates the allowed 
range from the global SM fit. 



fit to the four decay modes leads to the following results 



\T 

\P\ = 0.08 



0^53+0.029 



Q , |C| = 0.43 ± 0.05, Sc = -0.85 

+0.003 r _ n AQ+0.09 
-0.005' "P — 'J 



+0.52 
-0.28 

" -V - 1 11+"-" 
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(4) 



which shows a large ratio of \C/T\ ~ 0.81. In the 
QCD factorization estimation this value is bound to be 
|C/T| < 0.4. The error of Sc is significantly large. Note 
that the values of |T| and |P| and 7 remain almost un- 
changed, which indicates no explicit contradiction be- 



tween two sets of data tt^ 



~K and 7r°7r°, 7r°7r , and 



the relatively large ratio \C/T\ is purely the results of the 
large 7r°7r° and tt^tt" branching ratios. 

We now include other ttK data to determine Pew 
and its strong phase. We use the experimental value 
of Br{Tr°K°) = 11.5 ± 1.0, Br{-K- K^) = 24.1 ± 1.3 
and Br{TfK-) = 12.1 ± 0.8. The preHminary data of 
acp{Tr-K°) = 0.02±0.034 and acp(7r°if-) = 0.04±0.04 
are also considered. However, we do not include the pre- 
liminary data of acp(7r°7r°), acp{TT^ K'^) as we would like 
to leave them to be pure predictions from the fits. From 
the isospin structure of the effective weak Hamiltonian 
in the SM and the relations between the isospin ampli- 
tudes and the diagrammatic amplitudes, i.e., = Pew 
and 02 = T + C — Pew, one arrives at the following 
well-known model-independent constraint [lfij | 

Pew 3(6*9 + Cio) 



T + C 2{Ci + C2) 



(1.25 ±0.12) X IO- 



CS) 



with CiS being the Wilson coefficients evaluated at m?,. 
This relation tightly constrains the magnitude and the 
phase of Pew- However in the presence of new physics 
beyond the SM, the ratio could be significantly modified. 
In view of the recent puzzling experimental results, a 
careful analysis is urgently needed to find out whether 
this relation is indeed favored by the data. 



We now discuss several cases. First, consider a fit (Fit 
A) to the TTTT, ttK data using Eq.(|Sll. The result is given 
in the first column of TabQ] Comparing with the fit to 
7r+7r~ and 7r~^K~ in Eq.|21), one notices that the values 
of 7, |T| and \P\ are almost unchanged. C and its strong 
phase become larger and the ratio \C/T\ is enhanced to 
be close to ~ 0.9. Namely, the large tt^tt" branching ratio 
is actually not the full reason for a large C/T « 0(1). It 
is also required by the ttK data. The Xmin/'^-'^-f found 
to be 12.7/7 which is much higher than the previous fits. 
The main inconsistency comes from the branching ratio 
of TT^K~, 7r°A'° and 7r^/v°. The resulting best-fit val- 
ues in this case are 20.0 ± 0.8, 9.7 ± 0.5 and 22.3 ± 0.7 
respectively. The inconsistencies can be characterized 
by two ratios i?„ = Br{n+ K~)/ Br{Tr° K°) ~ 0.79 and 
R = Br{'n+K^)lBr{iT-K°) ~ 0.76, which should be 
very close to 1.0 in the SM. The small value of i?„ may 
require corrections to Pew while R may be connected 
to large non-facotrizable effects An important fea- 
ture of this fit is that the predicted direct CP violations 
acp(7r"7r") ~ 0.36 and acpiir^K^) ^ -0.11 _are large 
and compatible with acp(7r+7r^) and acpi'J^'^ K~). The 
Xmin 7 ^1^0 given in Fig^ which shows a good de- 
termination of 7. 
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FitB 


FitC 


FitD 
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1 0+"-" 

^•'-'-0.13 


J^-'J-o.is 


0.98+-" 


-1 1+0.12 

'-■'--0.19 


|T| 


0.52 ± 0.27 


0.52 ±0.03 


0.52 ±0.03 


1 10+0. 36 
^■^'-'-0.32 


\C\ 


0.47 ± 0.04 


0.45 ± 0.05 


0.45 ± 0.05 


"•'-'^-0.22 


Sc 


1 1+0.19 
-'-•-'--0.17 




1 07+0.3 


r, 7+1.29 
^- ' -0.3 


\P\ 


0.094 ± 0.001 


0.093 ± 0.002 


0.09 ± 0.002 


0.74 ±0.3 


Sp 


"-'•^-'-0.10 


-0.53iri:! 


-0.76 ±0.17 


n 0+0-05 

""•^-0.14 


\Pew\ 




0.03 ±0.01 


0.03 ±0.01 


0.024 ±0.0 






o.67i;- 


o.67i;;:i 


-,0+0.19 

^- ^^'-'-0.39 


\Po\ 


O(fix) 


O(fix) 


O(fix) 


0.11 ±0.02 


Spd 


O(fix) 


O(fix) 


O(fix) 


-0 21+" "^ 


ASp 


O(fix) 


O(fix) 


'^•^-0.17 


O(fix) 


xVdof 12.7/7 


9.1/5 


7.9/4 


5.4/3 




0.36 ±0.11 


0.05 ±0.22 


-0.06 ± 0.2 


0.07 ±0.39 




-0.10 ±0.004 


-0.01 ±0.05 


-0.02 ± 0.05 


-0.01 ±0.1] 




1.7 ±0.3 


1.56 ±0.4 


1.53 ±0.4 


1.7 ±0.5 




9.7 ±0.48 


11.1 ± 1.8 


11.1 ± 2.1 


11.3 ±2.3 




11.7±0.6 


11.9 ± 2.2 


11.8 ±2.4 


11.9 ±2.5 


«7r+7r- 


0.27 ±0.06 


0.30 ±0.08 


0.37 ±0.06 


0.34 ±0.27 




-0.1 ±0.02 


-0.11 ±0.02 


-0.1 ±0.03 


-0.1 ±0.06 



TABLE I: Best fitted parameters and predictions from charm- 
less B decay data in four different cases. Details are explained 
in the text. 

Second, considering a fit (Fit B) with freeing the pa- 
rameter Pew and its strong phase. The results are tabu- 
lated in the second column of tabled which show roughly 
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the same values of 7, |T|, \P\ and |C|, but the value of 
\Pew\ - 0.03 leads to 

(3.1 ± 1.3) X 10-2 (6) 

which is twice as large as in Eq.lQ. The data of 
Br{TT^K^) and the ratio i?„ are perfectly reproduced. 
This result agrees with the observation in Refs.0 with 
a statement that a large electro-weak penguin can con- 
sistently explain the ttK data. Clearly, the large value 
Pew is driven by the observed large branching ratio 
of tt^K" mode. All the previous fits with small Pew 
failed to meet this date point 0, Q. Note that in the 
case of large Pewi the predicted CP violations of tt'^tt" 
and -K^K^ are found to be small. The predicted central 
values are only 0.06 and —0.02 respectively though the 
errors are big(see TabP). Therefore it provides a possi- 
bility to distinguish the electro-weak penguin effects in 
the near future with more accurate measurements. In the 
third column (Fit C) of TabU we consider the effects of 
SU{?)) breaking in the strong phases. The best fit gives 
A(5p — 0.2 in accordance with Eq.Q. In this case, value 
of Pew remains the same as in Fit.B. The predictions 
give acp{T^^T^^) = —0.06 and acp(7r°A'°) = —0.1 respec- 
tively. 

The inclusion of all the tttt and nK modes allows one to 
investigate the possible large inelastic rescattering effect 
due to the process oi B —* DD(^s) 'K'k{ttK). It is well 
known that B DD have a large branching ratio about 
50 times greater than that of i? — > tttt, which ampli- 
fies the successive small effects of re-scattering I?-D(s) — ^ 
tt-kIttK). Considering the fact that B DD{s) only 
contributes to the isospin 0(1/2) irninK) final sates and 

is) 

carries only the CKM factor Ac , its contribution can be 
parameterized by only one complex quantity denoted by 
D{D') and effectively it can be considered by replacing 
pC) by P^^ = p(')-hD('). In the fourth column (Fit D) of 
TablU the parameters of |Pd| and 6p^ motivated by the 
inelastic rescattering from B —f DDi^g-^ are added which 
makes P and Pd two independent parameters. The re- 
sults show that Pd is compatible with the QCD penguins 
obtained from the previous fits of A,B,C while P becomes 
larger. This large difference between P and Pd indicates 
a large effect of inelastic rescattering and may also imply 
new physics in strong penguin sector. In this fit, the ratio 
of \C/T\ is reduced to 0.37. The ratio of Pew/{T + C) 
remains large and the two predicted CP violations are 
again small. 

In conclusion, the current data enable us to make a 
very encouraging global fitting for testing the standard 
model and probing new physics. It will be very crucial to 



arrive at more accurate measurements on both branch- 
ing ratios and direct CP violations in i? ^ tt^'tt^ and 
TT^K^. The current preliminary data give ac_p(7r"7r") = 
0.28±0.39 and acpiTT°K") = -O.O9±O.140. Due to the 
large errors, including them will not change the conclu- 
sion. Numerically, we find that the results in Eq.l^J are 
unchanged. The ratio \Pew/{T + C)\ remains large and 
is found to be 0.024 ±0.01, 0.034 ±0.01 and 0.033 ±0.04 
for FitB,C and D respectively in TabQI It is very likely 
that we are standing at the corner of finding new physics 
with two B-factorics. 
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